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ABSTRACT 

In an effort to explain the short-timescale variability of the broad, double- 
peaked profiles of some active galactic nuclei, we constructed stochastically per- 
turbed accretion disk models and calculated Ha line profile series as the bright 
spots rotate, shear and decay. We determined the dependence of the properties 
of the line profile variability on the spot properties. We compared the variability 
of the line profile from the models to the observed variability of the Ha line of 
Arp 102B and 3C 390.3. We find that spots need to be concentrated in the outer 
parts of the line emitting region to reproduce the observed variability properties 
for Arp 102B. This rules out spot production by star/disk collisions and favors 
a scenario where the radius of marginal self-gravity is within the line emitting 
region, creating a sharp increase in the radial spot distribution in the outer parts. 
In the case of 3C 390.3, all the families of models that we tested can reproduce 
the observed variability for a suitable choice of model parameters. 

Subject headings: accretion, accretion disks, line: profiles, galaxies: active 



1. Introduction 

1.1. Background and Motivation 

Double-peaked broad Balmer emission lines are found in 20% of radio loud active galac- 
tic nuclei (AGNs) at z < 0.4 (Eracleous & Halpern 1994; 2003 ) and 4% of the Sloan digital 



Sky Survey (SDSS) quasars at z < 0.33 (IStrateva et al. 20031 ). A number of models have 



been suggested for the origin of Double-peaked emission lines (DPEL), but most of them 
face significant challenges when compared to observations. The ve locity drift of the two 



peaks and the single-peaked profile of high-ionization emission lines ( lEracleous et al. 19971 ) 



are not consistent with the binary black hole scenario suggested by Gaskell (1983). The 
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observed response of the line profiles to changes in the ionizing flux (IDietrich et al. 1998 
OBrien et al. 19981 ) disagrees with predictions from the bipolar outflow model proposed by 
Zheng et al. (1990). The disagreement between the observed radio lobe geometry and the 
outflow geometry required to fit the line profile is also a challenge for this model. Wan- 
ders et al. (1995) suggested that the double-peaked profile could be produced by a col- 
lection of anisotropically illuminated clouds. Even though this scenario agrees with most 
observations, it has theoretical difficulties since these clouds would be destroyed rapidly 
by collisions or drag. Finally, Chen, Halpern, & Filippenko (1989) and Chen & Halpern 
(1989) suggested that the DPEL are emitted in a zone in an accretion disk over a nar- 
row range of radii (i.e. r out /r in < 10). In this scenario, the hi gh ionization emission lines 
are produced at the base of a disk- wind (jCollin-Souffrin 19871 ). Thus, the profiles of the 
high-ionization lines are single-peaked because of radiative transfer effects in the acceler- 
ating outflow ([Murray fc Chiang 19971 ). The limitation to this model is that the accretion 
disk structure required to explain the variability of the line profiles cannot be axisymmetric. 
However one does not expect a perfectly axisymmetric accretion disk so this is not a major 
challenge. The model of DPEL production in the accretion is the one that holds up best to 
close scrutiny, thus we adopt it as our working model. 

DPEL profiles are observed to vary on timescales of months to years, i.e. on timescales 
of the order of the dynamical time or longer (e.g., Veilleux & Zheng 1991; Zheng et al. 1991; 
Marziani et al. 1993; Romano et al 1998; Gilbert et al. 1999; Sergeev et al. 2000; Shapovalova 
et al. 2001; Storchi-Bergmann et al. 2003; Gezari et al. 2007). The line profile variability 
does not appear to be correlated to changes in the line and/or continuum flux so it likely 
traces changes in the accretion disk structure. The long-timescale variablity of the DPEL 



i a spiral arm ( 


Gezari et al. 2007: 


Gilbert et al. 1999 


, and references 



therein). These models, however, fail to explain t he long times cale variability of some objects 



and the short-timescale variability of all objects (ILewis 20051 ) 



Other attemps at explaining the line profile variability through perturbations of the disk 
structure have introduced bright spots over an axisymmetric accretion disk. Newman et al. 
(1997) succesfully modeled the variation of the Ha peak intensity ratio of Arp 102B with 
a single spot rotating within the disk, but Gezari et al. (2007) were not able to apply the 
same model to the same object at a different time period. It is noteworthy that a similar 
bright spot model was recently used by Turner et al. (2006) to explain the variability of 
the Fe Ka line profile of Mrk 766 in the X-ray band. Sergeev et al. (2000) tried a different 
approach; they modeled the variations of Arp 102B with a collection of 1500 clouds in 
a disk and were able to reproduce the root mean square (r.m.s.) spectrum and a quasi- 
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sinusoidal variation of the excess flux. This model however introduced phase effects and 
an edg e-on disk in disagreem e nt with the inclination obtained fro m fitting the total line 
profile (IChen fe Halpern 19891 ; IChen et al. 19891 ; iGezari et al. 20071 ). A collection of clouds 
orbiting the central black hole can also explain the power-law shape of the X-ray power 
spectra observed in AGNs (lAbramowicz et al. 1991c iBao fe 0stgaard 19951 ). Finally, Pariev 
& Bromley (1998) modeled the Fe Ka line profile resulting from a turbulent disk, which was 
created by adding a random frequency shift to each pixel in the accretion disk. Motivated 
by the ideas of Sergeev et al. (2000), we investigate in this paper whether a stochastically 
perturbed disk model can explain the short and long timescale variation of the DPEL Ha 
profiles of the two best-monitored objects: Arp 102B and 3C 390.3. 



1.2. The Physical origin of Bright Spots 

Many processes can lead to density or temperature inhomogeneities, hence brightness 
perturbations, in the accretion disk: self gravity, disk-star collisions, and baroclinic vortic- 
ity The properties of the inhomogeneities (size, radial distribution, lifetime, resistance to 
Keplerian shear) depend on their specific formation process, as we outline below. 

• A disk will fragment under the effect of self-gravity when the Jeans instability sets 
in. The radius of marginal self-gravity depends on the properties of the disk and the 
central object and can be as low as 500 r g (r g = GM/ c 2 is the gravitational r adius) 



for objects with an Eddington ratio of L/L-^dd ~ 10~ 3 ( jGoodman fe Tan 2004 ). The 
radius of marginal self-gravity can be even smaller with more extreme disk properties 
than those as sumed by these authors, su ch as a Shakura-Sunyaev viscosity parameter 



less than 0.3 (IShakura fe Sunyaev 19731 ). The size of the clumps will be smaller than 
the local Jeans radius rj = Cg/y/Gp oc £~ 3 / 20 , with a typical size of 10 r g at 1000 r g 
for a 10 8 M & black hole. The clumps produced in this manner do not shear with 
differential rotation and they have a high density (and hence brightness) that varies 
very little over time. The radial distribution of the clumps should be quite uniform 
past the radius of marginal self-gravity. 



Disk-star collisions are thought to occur as frequently as once per day (IZurek et al. 19951 ) 
Such collisions increase the local disk temperature, and hence create a bright spot. 
These spots will shear with Keplerian rotation and will fade away as the material cools 
down. The typical size of a spot at the time of formation is at least equal to the size 
of the star, which is very small in units of r g . As shown by Zurek et al. (1995), the 
radial distribution of the spots created by star collisions will increase with decreasing 
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radius since collisions are more frequent at smaller radii due to the increased Keplerian 
frequency of the star about the black hole: t / co iiison( r ) oc r p ~ 1 / 2 with < p < 1/2. 

The radial temperature gradient in the accretion disk combined wit h the Keplerian dif- 



ferential rotation of the disk material will lead to baroclinic vorticity (IPetersen et al. 20071 ) 



Gas drag will then cause the disk material to spiral to the center of the vortex increas - 



ing the local density, hence the brightness, of the vortex (IBarranco fc Marcus 20051 ). 
The lifetime and resistance of the spots to shearing is still a subject of debate. The 
simulations of Petersen et al. (2007) create long-lived, non-decaying, non-shearing 
vortices. Barranco & Marcus (2005) also find the same long-lived vortices when the 
vorticity is larger than the shearing; but they also find short-lived, shearing vortices 
when the vorticity is not high enough or when the gas speed is supersonic. Simulations 
by Johnson & Gammie (2005) lead to long-lived, non-shearing vortices that slowly 
decay as t~ 1 / 2 . All of these models produce structures whose size is comparable to the 
scale height of the disk, which scales as r 9//§ . The distribution of the spots produced 
by vorticity is weakly dependent on the radius. 

We investigate the effects of these different bright spot properties on the variability of 
DPEL profiles. In §2, we present a method to quantitatively characterize the observed line 
profile variability and apply it to two AGNs: Arp 102B and 3C 390.3. In §3, we describe 
the mathematical formalism for computing line profiles in our model. In §4, we explore the 
parameter space of the model and compare the results to the observations. In §5, we discuss 
the results and their implications . 



2. Charaterizing the Observed Line Profile Variability 

2.1. Available Data 

We use the spectroscopic observations of Arp 102B and 3C 390.3 described in Gezari 
et al. (2007). We chose these two broad-line radio galaxies because they were regularly 
monitored for extended periods of t ime and because their black hole masses are accurately 



measured (ILewis fc Eracleous 20061 ). The observations of Arp 102B were made from June 
1983 to September 2004, while the observations of 3C 390.3 span a shorter time period, 
from July 1988 to September 2004. These spectra have moderate resolution (4-6 A or 
180-260 km s" 1 ) and cover at least the wavelength range of the Ha line. 

After the initial processing of the data, the continuum and narrow emission lines were 
subtracted as described by Gezari et al. (2007). We excluded spectra for which Gezari et 
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al. did not achieve reliable subtraction of the narrow lines. We thus have 91 line profiles 
with an average time separation of 86 days for Arp 102B and 35 spectra every 174 days, on 
average, for 3C 390.3. Hence the dynamical timescale of the disk (as calculated in is well 
sampled. The average spectra of the two objects are shown in Figure [2] with the best fitting 
axisymmetric disk model superposed (see §1]). 

All these line profiles were normalized to a constant integrated flux. This was done be- 
cause the variability of the integrated emission line flux is correlated w ith the intensity of the 



ionizing continuum rather than with structural chan ges in the disk (e.g. iRosenblatt et al. 1992 



Kassebaum et al. 19971 ; IWanders & Peterson 19971 ). We then measure the seven line para- 



maters described by Gezari et al. (2007): the red and blue peak velocities (in km s _1 ), the 
blue peak to red peak flux ratio, the full width at half maximum and at quarter maximum 
(FWHM and FWQM, in km s _1 ) and the line displacement from its rest wavelength at 
half and quarter maximum (in km s _1 ) [see Figures 3 and 5 of Gezari et al. (2007) for the 
parameter time series of Arp 102B and 3C 390.3 respectively]. We chose to use these line 
parameters since they are independent of the normalization of the line profile and do not 
require any fitting of the line profile by a specific physical model. 



2.2. Method for Characterizing the Variability 

We cannot fit series of model line profiles to the observed data because of the inher- 
ently stochastic nature of the model (described in £J3]). Comparing the variation of the line 
profile parameters as a function of time is not meaningful either since a change in the initial 
observing time or initial conditions will modify the exact values of these parameters. Hence 
we developed a comparison method based on power spectra, which captures the statistical 
properties of the variability. 

Throughout this paper we use the Lomb-Scargle algorithm (Scargle 1982; Press & Ry- 
bicki 1989) to compute periodograms because our time series are unevenly sampled in time. 
We tested the robustness of the algorithm by comparing the periodograms of evenly- and 
unevenly-sampled, simulated time series. While the periodogram of an unevenly-sampled 
time series might have missing power peaks when compared to the results from an evenly- 
sampled series, it never has false positive power peaks. The sampling interval has a greater 
impact on the resulting power-law index of the periodogram: irregularly sampled time series 
produce significantly flatter periodograms because the Lomb-Scargle algorithm produces in- 
herently noisier periodograms at high frequency than the fast Fourier transform of an evenly 
sampled time series. To bypass this effect, we sample our simulated times series at the same 
times as the data. 
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From the line profile series, we created light curves in narrow velocity intervals (5X = 
20 A, corresponding to Av = 914 km s _1 ) and computed the periodograms. Thus, we 
produced 2-D periodograms, which are displayed in Figures [3] and HI for Arp 102B and 
3C 390.3 respectively. The probability that a power peak of height z is created at random 
is then P pa Ne~ z where iV is the number of frequency bins (110 and 36 for Arp 102B and 
3C390.3 respectively). Each 2-D periodogram has 40 velocity bins so one can expect the 
power to be above the 98% significance level in up to 88 bins for Arp 102B and 29 bins 
for 3C 390.3. The actual number of bins with power above the 98% significance level is 
59 for Arp 102B and 10 for 3C 390.3, which is below the random occurrence expectation. 
Moreover these power peaks are mostly located at low frequency where the error bars are 
large, making the significance of these peaks low. Because the significance of discrete power 
peaks is low, we concentrate on another aspect of the periodograms, namely the power-law 
index. The bottom panels of Figures [3] and H] show the collapsed (i.e., summed) periodograms 
as a function of frequency for all velocities v < (dashed line) and v > (continuous line). 
The power-laws that best fit each periodogram (assuming z oc f a ) are also displayed. For 
Arp 102B, the power-law of the periodogram with v < has an index of —0.60 ± 0.06 and 
that with v > 0, — 0.47±0.06. For 3C 390.3, the v < periodogram has a power-law index of 
—0.6 ±0.1 and the v > periodogram, —0.6 ±0.1. This indicates that the profile variability 
of Arp 102B is somewhat more asymmetric than that of 3C 390.3. The right panel shows 
the total power per unit frequency in each velocity interval. The velocity intervals with the 
highest average power are the most variable parts of the line profile. 

Another way to characterize the line profile variability is through the periodograms of 
the times series of profile parameters. From the line parameters that we previously mea- 
sured, we construct time series, and calculate the corresponding periodograms. The line 
profile parameter periodgrams for Arp 102B and 3C 390.3 are displayed in Figures [5] and [6] 
respectively. Power peaks seen in these periodograms are not considered significant due to 
their low power and their number which is consistent with chance occurrence. The power-law 
indices of all the line parameter periodograms of Arp 102B are negative and vary between 
—0.30 and —0.85. The line parameter periodograms of 3C 390.3 are flatter than those of 
Arp 102B, with power-law indices varying between —0.76 and 0. The values of all these in- 
dices are given in the captions of Figure [5] and [6J The error bars on these power-law indices 
are ~ 0.1, which means that: (a) the power-law indices of the periodograms of the red peak 
position and of the FWQM are significantly different from each other, (b) the power-law in- 
dices of the blue peak position periodogram is marginally different between the two objects, 
and (c) all the other power-law indices are consistent with each other. 

We also calculated the root mean square (rms) amplitudes of the line profile parameters 
for each object (listed in Tabled]). The fractional rms amplitudes (rms amplitude divided by 
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mean) of the line parameters give an alternative estimate of the amplitude of the variability. 
The fractional rms amplitude of the red and blue peaks positions and displacement at half 
and quarter maximum of 3C 390.3 are much larger than those of Arp 102B. This means that 
the line profile changes of 3C 390.3 have a much larger amplitude than Arp 102B. 

Hence the line profile variability properties of Arp 102B and 3C 390.3 are quite dif- 
ferent from one another whether we consider the 2-D periodograms, the line parameter 
periodograms or the fractional rms amplitude. With this in mind, we tested accretion disk 
structure models and their induced variability characteristics in search of a model that could 
reproduce the observed line profile variability of these two AGNs. 



3. Description of the Models 

We use the formalism laid out by Chen et al. (1989) and Chen & Halpern (1989) to 
calculate the double-peaked emission line profile from an accretion disk in the weak-field 
limit. We compute the line profile numerically by evaluating 



fu oc / U IuM, 4>, u) D 3 (£, </>) 0) d4 d£ 




over the surface of the disk. Here, <fi is the azimuthal angle in the disk frame and £ is 
the distance from the central black hole in units of r g . The velocity structure of the disk 
is contained in the Doppler factor D(£,<fi) and the gravitational redshift and light bending 
effects are included in \I/(£,0). These factors were computed by Chen et al. (1989) for a 
circular disk and depend on the inclination % of the axis of the disk to the line of sightfj]. 

The specific intensity I Ve (£, <f), v e ) contains information on the local line profile (assumed 
to be Gaussian) and the disk emissivity function e(£, 0): 



hMA,Ve) oc e(£,0)exp 



(^e - fo) 



21 



(2) 



2a 2 

where a is the line broadening parameter, u e is the emitted frequency of the photon and v Q 
is the rest frequency of the emitted photon. 

The stochastically perturbed model was implemented by creating the appropriate initial 
disk emissivity function e(£, <p). Over an underlying power-law emissivity of index — q (q > 0), 



1 The notation t/>(£, (f>) was introduced by Eracleous et al. (1995). This corresponds to the function g(D) 
used by Chen et al. (1989) 
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we add n circular spots. Each spot has a random position in the disk (£;, </>;), a radial Gaussian 
profile of spread Wi, and a contrast C-{. 

I »=i 

The values of & and are randomly generated. The radial distribution of the spots 
can be either uniform or a power-law, while their azimuthal distribution is always uniform. 
The spot size, Wi, scales as the disk scale height (oc £ 9//§ ) and the spot size at the inner 
radius (u> m in) is an adjustable parameter. The contrast, Ci, is randomly generated from a 
Gaussian distribution with a standard deviation of 1; the average value of Ci is an adjustable 
parameter. The average number of spots, n, is kept constant through a given simulation. 

We created three families of models from the basic spot distribution corresponding to 
the physical models discussed in §1. The properties of these families are as follows: 

1. The spots do not shear and do not decay as their centers rotate differentially around 
the central black hole. This model is relevant to disks in which self gravity is important 
or disks with long-lived non-shearing vortices, such as those described by Barranco & 
Marcus (2005) and Petersen et al. (2007). A snapshot of the disk as described by this 
model is shown in Figure UK 

2. The spots do not shear, but they do decay with time as their centers rotate differentially, 
which is relevant to the Johnson & Gammie (2005) model. A snapshot of the disk as 
described by this model is shown in Figure [Tb. 

3. The spots shear due to the Keplerian rotation of the disk and also decay. This model is 
relevant to spots created by disk-star collisions (power-law radial distribution of spots) 
and short-lived spots created by weak vorticity in the Barranco & Marcus (2004) model 
(with a uniform radial distribution of spots). A snapshot of the disk as described by 
this model is shown in Figure [lb. 

The Keplerian timescale for the rotation of the spots in the accretion disk depends on 
the mass of the black hole, Mbh, which can be a free parameter, in principle. The fading of 
the spots was treated by multiplying the initial value of C\ by an exponential decay function 
whose decay time Td eC ay is a free parameter. This decay time is assumed to be a function of 
radius (oc £ 3 / 2 ), i.e. proportional to the local dynamical time. When the spots decay away, 
new spots have to be created to keep the average number of spots constant. This is done by 
balancing the average decay rate and the appearance rate. 



2w? 





(e - &r 


exp 


2wf 





(3) 
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Hence the free parameters are: q, M B h, £min, U> h (Ci)> ^min, n, and r decay , if 
the spots decay. The black hole masses are set equal to the values reported by Lewis & 
Eracleous (2006) and we assume that q = 3, following Dumont & Collin-Souffrin (1990). 
Over the total baseline spanned by the observations, the inner and outer radius of the 
line-emitting portion of the accretion disk could vary. For example, the inner radius could 
increase at a time of increased l uminosity of the central engine due to ionization of the 



inner parts of the accretion d isk (jCollin-Souffrin fc Dumont 19901 ; IShapovalova et al. 2001 



Storchi-Bergmann et al. 20031 ). If we let £ min and £ max be free parameters for each profile 
of a time series, the total number of free parameters increases greatly (e.g., 185 or 186 
free parameters for Arp 102B). Instead we use an average value of £ min and £ max obtained 
from fitting an axisymmetric model to the average and minimum spectra. The minimum 
spectrum was constructed by taking the minimum value of the normalized flux from the 
entire time series at each wavelength. Keeping the radii fixed is justified by the fact that 
large-amplitude continuum variability (changes by more than an order of magnitude) occurs 
on timescales considerably longer than the Keplerian timescale at these radii. Therefore 
noticeable fluctuations in the radii should also occur on timescales considerably longer than 
the Keplerian timescale at these radii. This reduces the total number of free parameters to 
3 or 4. By keeping the radii fixed we are able to isolate the effect of emissivity perturbations 
on the variability of the line profiles. 

The computational implementation of the model requires dividing the disk into a loga- 
rithmic grid in £ and linear grid in 0. We used a 200 x 200 grid in order to resolve each spot 
into more than a few pixels and to resolve the effect of Keplerian shear. We chose a logarith- 
mic division of the grid in the radial direction because the bulk of the line flux is produced 
at small radii (due to the £ -3 emissivity function) and because most of the quantities that 
depend on radius are described by power laws. 



4. Application 
4.1. Arp 102B 

We constructed the minimum and average spectra from the 91 line profiles of Arp 102B 
and fitted both of them with a model employing an axisymmetric power-law emissivity of 
the form e oc with q = 3 and i = 30°. The best fit results for the average profile are 
Um = 290 ± 30, £ max = 1000t|$ and a = 1300 ± 400 km s _1 . The best fit parameters to 
the minimum spectrum are £ min = 340 ± 40, £ max = 940^2oo an d cr = 1200 ± 400 km s _1 , 
which are quite consistent wit h the results of the best fit to the averag e spectrum and wit h 
fit results by previous studies (jchen fc Halpern 1989 ; Chen et al. 1989 ; Gezari et al. 2007 ). 
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The average spectrum and the best fit are shown in Figure^. 

For the analysis of Arp 102B we set the underlying disk parameters to the following 
values: £ min = 300, £ max = 950 and a = 1200 km s _1 . The mass of the black hole is set to 
10 8 M Q as measured by Lewis & Eracleous (2006). Hence the dynamical times (Td yn ) at the 
inner and outer disk radii are 0.08 and 0.45 years, respectively. 

We note that the fit to the average spectrum is not as good as what can be obtained 



for individual spectra (e.g., IChen fc Halpern 19891 ) . We investigated the effect of setting 



the disk parameters to different values on the resulting variability properties. We used the 
best fit values found by Chen &; Halpern (1989), namely i = 32°, £ min = 350, £ out = 1000 
and a = 850 km s -1 and the values we decided upon in two runs of the first family of 
models with identical spot properties. We then computed the line parameter periodograms 
and found that the difference in power-law indices between the two runs was very small 
Aa = 0.08 ± 0.04. This is a result of the subtraction of the average profile before computing 
the periodogram. Since we are concerned with perturbations about the mean, the exact 
values of the disk parameters have a negligible effect on our conclusions. 



4-1.1. First Family of Models: Non- Shearing, Non-Decaying Spots 

We first characterize the variability of the model with non-shearing, non-decaying spots 
(i.e. the first family of models). In order to investigate the effects of varying the free 
parameters, we vary only one parameter at a time (namely, the number, size, contrast, and 
radial distribution). We then create 100 time series with randomly selected values of the 
parameter of interest, measure the power-law indices of the line parameter periodograms 
for each run, and plot the average and standard deviation of the power-law indices against 
the free parameter value. Figure [7] shows as an example the dependence of the power-law 
indices of the line parameter periodograms on the number of spots with the other spot 
parameters fixed to w min = 10 r g , (Cj) = 10 and a constant number of spots per unit 
radius. We calculated the Pearson coefficient and the significance of the Spearman rank 
coefficient in order to search for any correlation between the parameter of interest and the 
power-law indices. Since the dispersion of the power-law indices is large, we set stringent 
limits for considering these coefficients to be significant (||Pe|| > 0.8 and P s < 0.02). We 
also examined the plots to look for a break in the correlation and recalculate the coefficients 
on both sides of the break when we suspect one. Varying the number of spots produces 
power-law indices consistent with or greater, which does not match the observed values 
(—0.30 and —0.85). There is no trend between any power-law index and the number of spots 
in this model. Using the same approach for the other model parameters of interest, we find 
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that the logarithmic slope of the power spectra depends only weakly, if at all, on the model 
parameters. The weak trends that we notice are: 

• All the periodograms become steeper as the spot contrast increases. 

• Most of the periodograms become steeper when the spots are concentrated in the outer 
parts of the disk. 

• The power-law index of the periodogram of line shifts at QM has a broad maximum 
at a spot size ~ 10 r g ; it is less than 0.3 for small (< 5 r g ) and large (> 12 r g ) spots 
while it is ~ 0.4 around 10 r g . 

Thus bright spots with W7 min < 5 r g or w m i n > 12 r g , concentrated in the outer parts 
of the disk produce the steepest possible power spectra. However, even the steepest model 
power spectra are too flat to explain the power spectra observed in Arp 102B (see illustration 
in Fig. ED. 

We created 2-D periodograms for each run of the model with non-shearing, non-decaying 
spots and calculated the power-law indices of the collapsed periodograms for v < and v > 0. 
The difference between the two power-law indices vary from —0.2 to 0.2. The average value 
for the power-law indices is —0.2, which is higher than the observed values for Arp 102B. 
The power-law indices have a slight dependence on the radial spot distribution: the most 
negative power-law indices are produced when spots are concentrated in the inner part of 
the disk. 

We also investigated the effects of varying the model parameters on the fractional rms 
amplitude. We find that the fractional rms amplitude increases with increasing contrast and 
size of the spots, but does not depend significantly on the number of spots or the radial 
distribution. Furthermore, in order to reproduce the observed fractional rms amplitude of 
Arp 102B, a combination of particularly large size and contrast is necessary (C; > 20 and 
w min > 10 r g ). 

In conclusion, the first family of models can produce asymmetric 2-D periodograms, as 
observed. However, this family also produces power-law indices that are too high compared to 
what is observed in Arp 102B (i.e. the simulated power spectra are too flat). The agreement 
between this family of models and the observed properties of the line profile parameters of 
Arp 102B is best when the spots are concentrated in the outer parts of the disk. But the 
best agreement with the collapsed periodograms is obtained when spots are concentrated in 
the inner parts of the disk. 
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4-1-2. Second Family of Models: Shearing, Non-Decaying Spots 

Similarly, we characterize the variability of models whose spots decay and do not shear 
(i.e., the second family of models) by letting only one parameter free to vary and by measuring 
the spread in power-law indices in each parameter bin. We find that: 

• The periodograms of the line shifts at QM and HM and of the red peak position steepen 
when the decay time gets shorter. However the power-law indices of the FWHM and 
FWQM steepen when the cooling time gets longer. 

• The periodograms of the blue peak position steepen with decreasing initial contrast of 
the spots. 

• The periodograms steepen when the radial distribution of the spots favors more spots 
in the outer part of the disk. This correlation is not as strong as for the model without 
spot decay (first family of models). 

Compared to the first family of models, the second family produces periodograms whose 
power-law indices are steeper by 0.2 for a decay time of 5rd yn and 0.4 when Td ccay = Td yn . This 
brings the simulated power-law indices closer to the observed values. The power-law indices 
of the periodograms of the peak flux ratio, blue peak position, the FWHM and the FWQM 
are in good agreement with the observed values of Arp 102B for a £ 2 radial distribution of 
spots and r dccay > 5r dyn , but the power-law indices for the other line profile parameters are 
too high by 0.3 to 0.6 

The difference between the two power-law indices of the collapsed line flux periodograms 
varies between —0.2 and 0.3. The average value for the power-law indices is 0, which is greater 
than for the first family of models and the observed values. The power-law indices show the 
same dependence on radial spot distribution. The fractional rms amplitude is similar to that 
of the first family and the decay time has no effect on the fractional rms amplitude. 

4-1.3. Third Family of Models: Shearing, Decaying Spots 
Finally, for the model whose spots shear and decay (i.e. the third family of models), we 

find: 

• The periodograms of the FWHM and FWQM steepen with increasing decay time, 
increasing contrast, and decreasing number of spots. 
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• The periodograms of the red and blue peak position steepen with decreasing contrast. 

• The periodograms steepen when the radial distribution of the spots favors more spots 
in the outer part of the disk. The strength of this correlation is comparable to that 
observed for the second family of models. 

This family of models produces periodograms with power-law indices that are exactly 
like those of the second family of models, within the dispersion. The 2-D periodograms 
produced by this family of models have the same behavior as those of the second family of 
models. The fractional rms amplitude is similar to that of the first and second families of 
models and the decay time has no effect on the fractional rms amplitude. 



4.2. 3C 390.3 

Following the same procedure as for Arp 102B, we constructed the mean and average 
spectra from the 35 line profiles of 3C 390.3 and fitted both of these spectra with an axisym- 
metric power-law emissivity model of the form e oc £ -3 . The best fit results for the average 
profile are i = 27° ± 2°, £ min = 420+gJ, £ out = 1400^, a = 1500 ± 300 km s _1 . The best 
fit parameters to the minimum spectrum are i = 27° ± 1°, £ m i n = 500^4q°, £ out = 1500+3qq°, 
a = 1200^4qq km s™ 1 . The average spectrum and the best fit are shown in Figure [2b. 

Thus, for the analysis of 3C 390.3 we set the disk parameters to the following values: 
i = 27°, £ m i n = 450, £ out = 1400 and a = 1300 km s _1 . The mass of the black hole is set 
to 3 x 1O 8 M (see Lewis & Eracleous 2006). Hence, the dynamical times at the inner and 
outer radii are 0.45 and 2.5 years, respectively. 

We explored the parameter space of the three families of models as we did for Arp 102B 
( §4.ip and we found similar correlations for all the families. The correlations of the power-law 
indices with the radial distribution of spots are not as strong as those obtained for Arp 102B. 
This is probably a result of the fact that 3C 390.3 has a larger disk size and naturally has 
more spots at large radii than Arp 102B. 

Since the observed power-law indices of the periodohrams of 3C 390.3 are lower than 
those of Arp 102B, a wide range of model parameters can reproduce the values observed for 
3C 390.3. Figure [5] shows that the power-law indices of simulated periodograms are close to 
the observed values for 3C 390.3, here for a shearing spot model with a £ 2 radial distribution 
of the spots and variable r^ ecay . Such a close match is easily obtained for any family of 
models, depending on the initial parameter values. 

The 2-D periodograms of 3C 390.3 have the same behavior and same typical values of the 
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power-law indices as those of Arp 102B. Similarly, the fractional rms amplitude increases with 
increasing contrast and size of the spots, but is not sensitive to the other model parameters. 
A combination of large spot contrast and large spot size is required to match the observed 
fractional rms amplitide of 3C 390.3. 

5. Discussion 

Stochastically perturbed disk models produce variability of the Ha line profile whose 
properties depend on the model parameters. Allowing for decay of the spots steepens the 
observed periodograms of the line parameters. This is because the decay introduces a new, 
longer timescale in the system, thus shifting power to lower frequencies The radial distribu- 
tion of the spots also influences the power-law indices of the line parameter periodograms: 
steeper periodograms are produced when spots are concentrated towards the outer parts of 
the line-emitting region. This is a result of the fact that in such a distribution more spots 
have longer orbital periods, introducing more power at lower frequencies, thus steepening 
the periodogram. Shearing does not significantly change the power-law indices of the peri- 
odograms because shearing and decaying happen simultaneously and the effect of decaying 
of a spot is more prominent that that of shearing. Because shearing does not modify the 
variability properties of the line profile, one cannot constrainr it from observations of the 
line profile variability. The fractional rms amplitude of the line parameters increases with 
increasing contrast and size of the spots and does not depend strongly on the number of 
spots, their radial distribution, or their shear and decay properties. The fractional rms am- 
plitude thus provides another constraint on the contrast of the perturbation of the emissivity 
function. 

The perturbed disk model can reproduce the observed variability properties of Arp 102B 
and 3C 390.3 with different levels of success. Most of the observed line parameter peri- 
odograms of Arp 102B are too steep to be successfully reproduced. The best models have 
decaying spots, mostly concentrated towards the outer parts of the line emitting region. Such 
models produce power-law indices of the peak flux ratio and blue peak position periodograms 
that are comparable to those of Arp 102B, but all the other line parameter power-law indices 
are too high by up to 0.4. The maximum radial spot distribution index that we tested is 
+2. A model with an even steeper radial distribution might produce steep enough peri- 
odograms to match the observations. Such a distribution rules out the production of spots 
by collisions with stars since this process skews the spot distribution towards the inner part 
of the accretion disk. Spot production by vortices yields a priori a uniform spot distribution. 
Spot production by self gravity leads to spots only past the self gravity radius so this could 
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skew the spot distribution towards the outer part of the line emitting region. These spots 
would however neither shear nor decay so one would need an even steeper radial distribution 
index to achieve agreement with the observed values of Arp 102B. A sharp increase from 
to a large number of spots past the radius of marginal self-gravity might mimic such a 
steep radial distribution index. Using the properties of Arp 102B determined by Lewis & 
Eracleous(2006) and the formulae for the radius of marginal self-gravity (r sg ) given by Good- 
man k Tan (2004) and Hure (1998), we estimate that 400 < r sg /r g < 900 by varying the 
disk parameters within their plausible range. This estimate falls within the range of the line 
emitting region determined by fitting the average spectrum. To test whether a sharp bound- 
ary in the spot distribution steepens the power-law indices of the periodograms, we ran one 
simulation of the first family of models with no spots below 750 r g and 100 spots uniformly 
distributed above that radius. The spots were large (15 r g ) and bright ((Cj) = 20) as one 
would expect from self-gravitating clumps. The resulting line parameter periodograms had 
power-law indices varying between —0.17 and —0.72, which is steeper than what a smooth 
radial spot distribution produces and closer to the values observed for Arp 102B. Moreover, 
these self-gravitating spots would have the necessary size and contrast to match the observed 
fractional rms amplitude of Arp 102B. This result bolsters the hypothesis that spots inhabit 
the outer parts of the line-emitting disk of Arp 102B. 

There could also be other factors that steepen the line parameter periodograms. We 
assumed constant inner and outer radii of the line emitting region, but we tested the effect 
of a fast modulation of these radiii and found that this produces steeper periodograms. The 
inner and outer radii could vary due a change in the disk ionization structure after a variation 
of the ionizing continuum from the central source. Developing a physically self-consistent 
model to account for this effect requires one to calculate the variation of the whole disk 
structure with ionizing flux and is beyond the scope of this paper. 

The line parameter periodograms of 3C 390.3 are not as steep as those of Arp 120B 
and can be repoduced by a wide range of model parameters. A radial distribution of spots 
skewed towards the outer parts of the line emitting region is not absolutely necessary for 
3C 390.3 for two reasons. First, the steepening of the line parameter periodograms with 
increasing radial spot distribution index is not as marked as for Arp 102B. This might be a 
result of the fact that the line emitting region of 3C 390.3 has a much larger outer radius 
than that of Arp 102B. Hence, there is naturally a larger number of spots at large radii. This 
suggests that it is the relative number of spots at large radii compared to lower radii rather 
than the absolute number of spots at large radii that influences the power-law indices of the 
periodograms. Also since the observed power-law indices of the line parameter periodograms 
of 3C 390.3 are higher than those of Arp 102B, we do not need to select the model families and 
parameters producing steep power spectra in order to reproduce the properties of 3C 390.3. 
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This also means that no particular model for the origin of spots in the accretion disk of 
3C 390.3 can be ruled out or favored. We note that 250 < r sg /r g < 400 for 3C 390.3, which 
is below £ min and strengthens our assumption of a perturbed accretion disk at all radii. 

The behavior of the collapsed periodograms from the 2-D periodogram is a greater 
challenge to this model. The difference between the two power-law indices agree with the 
observed difference, but the values of the indices are too high to be reconciled with the 
observations. 

The results of our line profile variability study are quite different between Arp 102B 
and 3C 390.3 even though their black hole masses are comparable. The main driver of this 
difference is the higher accretion rate of 3C 390.3 (0.15 M Q yr _1 vs. 1.7 x 10~ 3 M yr _1 for 
Arp 102B), which makes that disk more unstable to self-gravity. It is encouraging that the 
perturbed accretion disk model provides a unified solution to explain the profile variability 
of two AGNs that are so different. In the context of the production of clumps by self-gravity, 
AGNs with similar black hole masses and even higher Eddington ratios than 3C 390.3 are 
expected to be unstable at all radii and the variability of their profile would then be well 
explained by our model. AGNs with smaller black hole masses might have their radius of 
marginal self-gravity in the line-emitting region and have a variability pattern similar to 
Arp 102B. In order to test this theory and confirm our result, we need to apply this method 
to a larger sample of AGNs covering a wider range of black hole masses and accretion rates, 
such as those presented in Gezari et al. (2007) and Lewis (2005). 

This work was partially supported by the Zaccheus Daniel grant. We would like to 
thank Suvi Gezari for providing fully reduced spectra for Arp 102B and 3C 390.3. We also 
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Table 1: Observed fractional rms amplitude of the line profile parameters of Arp 102B and 
3C 390.3 







Arp 102B 


3C 390.3 


Red Peak Position (km s" 




(8.9 ±0.3) x 10" 2 


0.35 ± 0.03 


Blue Peak Position (km s" 


- 1 ) 


(4.4 ±0.4) x 10~ 2 


(6.7 ±0.1) x 10~ 2 


Peak Flux ratio 




0.10 ±0.02 


0.15 ± 0.02 


FWHM (km s- 1 ) 




(3.9 ±0.2) x 10~ 2 


(3.8 ±0.03) x 10~ 2 


FWQM (km s- 1 ) 




(4.4 ±0.2) x 10~ 2 


(4.3 ±0.02) x 10~ 2 


Displacement at HM (km 




0.40 ± 0.04 


4.6 ± 0.6 


Displacement at QM (km 


s- 1 ) 


0.31 ± 0.03 


1.0 ± 0.1 
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Fig. 1. — Snapshots of the disk emissivity function for the stochastically perturbed disk 
model without decay or shearing of the spots (top), with decay but without shearing (middle), 
and with decay and shearing (bottom). Each simulation was started with identical initial 
conditions and the snapshots were all made 6 months after the initial time. Darker shades 
represent areas of higher emissivity. 
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Fig. 2. — Average spectra of Arp 102B and 3C 390.3 with the best fitting axisymetric disk 
model superposed as a dashed line. The model is described in §3] and the values of the model 
parameters are given in sections §4.11 and 14.21 
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Fig. 3. — 2D periodogram of Arp 102B computed in 20 A bins. The grayscale indicates 
the power. The horizontal dashed line represents the rest wavelength of the emission line. 
The vertical dashed line represents the Keplerian frequency of the outer radius of the line 
emitting region. The total power per unit frequency in each velocity bin is displayed in the 
right panel. The bottom panel shows the total power per unit frequency as a function of 
frequency, summed over all bins with v > (continuous) and v < (dashed), along with the 
best fitting power-laws (assuming z oc f a ). The indices of these power-laws are — 0.47±0.06 
and —0.6 ±0.06, respectively. (Note that the frequency has a logarithmic scale in the bottom 
panel and a linear scale in the top panel.) 
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Fig. 4. — Same as Figure [31 but for 3C 390.3. The left and right vertical dashed lines in the 
top panel represent the Keplerian frequencies of the outer and inner radii of the line-emitting 
part of the disk, respectively. The power-law indices are —0.6 ±0.1 and —0.6 ±0.1 for v < 
and v > 0, respectively. (Note that the frequency has a logarithmic scale in the bottom 
panel and a linear scale in the top panel.) 
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Fig. 5. — Periodogram of the line profile parameters of Arp 102B. The horizontal dash-dot 
line indicates the 98% signifance limit of the power peaks. The slanted dashed line shows 
the best-fitting power-law (assuming z oc f a ). The power-law indices are (left to right, then 
top to bottom) -0.50, -0.30, -0.45, -0.80, -0.70, -0.85, -0.85 with an uncertainty of 
~ 0.1. The vertical dashed line indicates the Keplerian frequency of the outer radius of the 
line-emitting part of the disk, obtained from the model fits described in §H 
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Fig. 6. — Same as Figure El but for 3C 390.3. The left and right vertical dashed lines 
represent the Keplerian frequencies of the outer and inner radii of the line-emitting part of 
the disk, respectively. The power-law indices are (left to right, then top to bottom) —0.07, 
-0.49, -0.36, 0.00, -0.59, -0.73 and -0.76 with an uncertainty of ~ 0.1. 
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Fig. 7. — Average and standard deviation of the power-law indices of the line parameter 
periodograms as a function of the number of spots in the disk. The vertical lines represent 
the scatter in the power-law indices of the simulations, not the error bars of the measured 
power-law indices. The model used here includes non-decaying, non-shearing spots. We 
investigate the effect of varying the number of spots between 50 and 500, while keeping all 
other parameters constant. The horizontal dashed lines represent the value of the power- law 
index found in the corresponding observed periodogram of Arp 102B. 
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Fig. 8. — Examples of realization of the shearing spot model for which the power-law indices 
of the line parameter periodograms are close to the observed values for 3C 390.3. The ratio 
of the decay time to dynamical time is the only parameter that was allowed to vary. The 
horizontal dashed lines represent the value of the power-law index found in the corresponding 
observed periodogram of 3C 390.3. 



